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, Abstract 



We study rates of convergence of solutions in L 2 and H 1 / 2 for a family of elliptic 
systems {C £ } with rapidly oscillating coefficients in Lipschitz domains with Dirichlet 
or Neumann boundary conditions. As a consequence, we obtain convergence rates 
for Dirichlet, Neumann, and Steklov eigenvalues of {C £ }. Most of our results, which 
rely on the recently established uniform estimates for the L 2 Dirichlet and Neumann 
problems in [12\ I13j. are new even for smooth domains. 

1 Introduction 

Let u £ G H 1 (Q) be the weak solution of C £ (u £ ) = F in Q subject to the Dirichlet condition 
u £ = f on dtt, where F G L 2 (tt), f G # 1/2 (<9ft) and C £ = -div[A(x/e)V]. Assuming that 
the coefficient matrix A(y) is elliptic and periodic, it is well known that u £ — > uq weakly in 



H 1 ^) and strongly in L 2 (Q), where uq G is the weak solution of the homogenized 

system Cq{uq) = F in Q and uq = f on dfl (see e.g. [1]). The same holds under the Neumann 
boundary conditions |^ = |^ = g G H^ 1 ^ 2 (d^l) with < g, 1 >= — J n F, if one also requires 
In Ue = In u ° = ®- ^he P r i mar y purpose of this paper is to study the rate of convergence 
of \\u £ — u \\ L 2(ty, as e — > 0, in a bounded Lipschitz domain Q C M d . As a consequence, 
we obtain convergence rates for Dirichlet, Neumann, and Steklov eigenvalues of C £ . Most 
of our results, which rely on the recently established uniform regularity estimates for the L 2 
Dirichlet and Neumann problems in (T2J HB) . are new even for smooth domains. 
More precisely, we consider a family of elliptic systems in divergence form, 

e>0. (1.1) 



c- 9 



dxi 



a p fx\ _d_ 
a v \e) dx< 



We will assume that A{y) = (a"-(y)), 1 < i, j < d, 1 < a, (3 < m is real and satisfies the 
ellipticity condition, 

m 2 < <${V)&$ < -Kl 2 for y G R d and £ = (£f) e M dm , (1.2) 
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where /x > 0, and the periodicity condition 

A(y + z) = A(y) for y E R d and z E Z d . (1.3) 

We shall also impose the smoothness condition, 

\A(x) - A(y)\ < t\x - y\ x for some A e (0, 1) and r > 0, (1.4) 

and the symmetry condition A = A*, i.e., aff(y) = aj°*(y) for 1 < i, j < d and 1 < a, {3 < m. 
We say A E A(/x, A, r) if it satisfies conditions ( fT2|) . ( TOj) and ( Oil . 
The following are the main results of the paper. 

Theorem 1.1. (Dirichlet condition) Let Q be a bounded Lipschitz domain, A E A(/x, A, r) 
and A* = A. Gzven F e L 2 (fi) and f E H 1 (d^l), let u £ G ff'(fl), £ > be the unique 
weak solution of the Dirichlet problem: C £ (u £ ) = F in Q and u e = f on dQ. Then for 
0<e< (1/2), ' 

IK - wolk 2 ^) + \\M(u e - m )||l 2 (oq) < C e\\u \\ H 2 {n) (1.5) 

if u E H 2 (Vt), and 

\\u e -u \\mn) < C a e\\ia(e)\^ + ' 7 {\\F\\ L 2 (n) + \\f\\ H Hdn)} , g . 
||.M(w £ - wo) ll^an) < CU| ln( £ )|i +CT {||F|| L2(Q) + \\f\\m(an)} 

for any a > 0. 

Theorem 1.2. (Neumann condition) Let Q be a bounded Lipschitz domain, A E A(/x, A, r) 
and A* = A. Given F E L 2 (Vt) and g E L 2 (dVt) with J n F + J m g = 0, letu £ E H 1 ^), e>0 
be the unique weak solution of the Neumann problem: C £ (u £ ) — F m O, ^f- = g on d£l and 
f n u £ = 0. Then for < e < (1/2), estimate U.5\) holds if Wo E H 2 (Q), and 

IK - w Q || L 2(n) + ||w £ - wo||L2(9fi) < C CT ^| ln(e) | ^ +<T {||F|| L 2 (n) + ||#||L2(an)} , , . 

3 \ / 

\\M(u e - w ) |U2(an) < C CT e| ln(e)|5 +<T {||F|| £ 2 (n) + y^en)} 

/or any a > 0. 

Here and thereafter M. denotes the radial maximal operator (see Section 2 for its defini- 
tion). Note that estimates of M.(u £ — Uq) in L 2 (dQ) in Theorems II. llTOl imply, in particular, 
the convergence of u £ to Uq in L 2 (S) uniformly for any "parallel boundary" S of Q. Also 
observe that in the case F = 0, estimates (11.6j) - (jl.7j) give 

~ u \\ L 2(n) < C a e\ \n(e)\^ +a \\uo\\ H i {dQ) , (1.8) 

for any cr > 0. The estimate of u £ — uq in L 2 (Q) by «o and its first-order derivatives is a 
natural question in the theory of homogenization (see [8] for a two-dimensional result). It is 
not known whether the logarithmic factor in fjl.Sp is necessary, even for smooth domains. 

We now describe the existing results on L 2 convergence and our approach to Theorems 
II. 1111.21 For a single equation (m = 1) with the Dirichlet condition, it is known that 

IK - w ||L2 (n) < Ce {\\V 2 u \\ L 2( Q) + ||Vn || L oc ( ^)} , (1.9) 
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holds without any smoothness or symmetry condition on A(y) or smoothness of Q. To see 
this, one considers 



w e (x) = u e (x) — Uq(x) — ex{x/e)VuQ{x) in Q, (1.10) 

where x{v) ls the matrix of correctors for C £ . Let w e (x) = 6 e (x) + z £ (x), where 9 £ is 
the solution to the Dirichlet problem: C £ {9 £ ) = in Q and 9 £ = — ex{x/e)Vu on dVL. It 
follows from the energy estimates that || Vz e ||#i(Q) < Ce\\ V 2 uo||,L 2 (r2) (see e.g. [T0l[T5]). This, 
together with the estimate ||# e ||L°°(Q) < Ce\\ VMo||L°°(sn) obtained by the maximum principle, 
gives f ll.9p . For elliptic equations and systems in a C 1,a domain with A G A(/x, A, r), the 
uniform estimates in [2J, [TJ for the L 2 Dirichlet problem imply ||0 E ||.L2(n) < Ce\\ ^7uo\\L 2 (dn) < 
Ce||?xo|| j? 2 (s^)- It follows that 

\\u £ - < H^Hl 2 ^) + ||#e|U 2 (^) + C^||VMo||L2( n ) < Ce||Mo||fl-2(J2), 

as noted in [15J. Using the recently established uniform L 2 estimates in [13], in the presence 
of symmetry (A = A*), we extend this result to the case of Lipschitz domains in Section 3, 
where we in fact prove that 

\\M(w £ )\\ L 2 (m) + \\w £ \\ H i/2 (n) + jy \Vw £ (x)\ 2 5(x) cfej < Ce||uo||jffa(n), (1-H) 

where 5(x) = dist(x, dfl) (see Theorem 13.41) . and deduce (11.51) as a simple corollary of (11. lip . 

The proof of (I1.6P is more involved than that of (11.51) . Note that with boundary data 
/ G iJ 1 ((9r2), one cannot expect Uq G H 2 (Q). Furthermore, if Q is Lipschitz, uq may 
not be in H 2 (Q) even if F and / are smooth (it is known that Uq G H 3 ^ 2 (Q) |9J). To 
circumvent this difficulty, our basic idea is to replace uq in (ll.lOp by a solution v £ to the 
Dirichlet problem for Co in a slightly larger domain: Cq(v £ ) = F in Q £ and v £ = f £ on 
dQ £ , where Q £ is a Lipschitz domain such that Q £ D Q and dist(<9fi e , <9fi) ~ e. Also, F 
an extension of F and f e (Q) = f(Aj 1 (Q)), where A e : dfl —> dQ £ is bi-Lipschitz map. Let 
w £ = u £ — v £ — ex(x/e)Vv £ = z £ + 9 £ , where 9 £ solves 

C £ {9 £ ) = in n, , 
J y ' (1.12) 
9 £ = f — v £ — ex{x/e)Vv £ on dVL. 

The desired estimates of 9 £ follow from the estimates for the L 2 Dirichlet problem in [13J. To 
handle z £ , one observes that £ £ (z £ ) = e div(h e ) in Q and z £ = in dQ, where \h £ \ < C|V 2 t> £ | 
in Q. Using weighted norm inequalities for singular integrals, we are able to bound ||£ e ||i,2(Q) 
and ||A^(^)|| L a(an) as wel1 as ll%llsva(n) h Y 

Ce I \V 2 v £ (x)\ 2 5(x)MS(x)) rfx} ' < Ce\ m(e)|f {\\F\\ L2{Q) + \\f\\ HHa a)} , 

for suitable choices of a's, where (j) a (t) = {ln(~ + e a )} a . See Section 4 for details. 

Very few results are known for the convergence rates in the case of the Neumann boundary 
conditions. By multiplying by a cut-off function the third term in the right hand side of 
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fll.lOp . one may obtain an 0(y/e) estimate of \\u E — Wo||L 2 (n), regardless of the boundary 
condition [U HP] . As far as we know, the only other known result is contained in [16], where 
the estimate \\u £ — Mollis) < C£||wo||h2(q) was proved in a curvilinear convex polygon Q 
in M 2 . In Section 5 we prove estimate (jl.5p in bounded Lipschitz domains in M. d , d > 2 for 
the Neumann boundary conditions. The proof uses an explicit computation of the conormal 
derivative ^f- on dVl and relies on the uniform estimates for the L 2 Neumann problem in 
[T2"j [T5] . The proof of estimate (jl.7p . which is given in Section 6 and also uses estimates for 
the L 2 Neumann problem in [TJl US], is similar to that of (jl.6p . It is interesting to point out 
that in this case the function v E , which replaces Uq in (11.101) . is a solution to the Dirichlet 
problem for £ m ^e, with boundary data given by a push- forward of uo\an- 

By a spectral theorem found in [10J , the L 2 error estimates of u £ — uq in Theorems 11.1111.21 
lead to error estimates for eigenvalues of {C £ }. For e > 0, let {fi £ } denote the sequence of 
Neumann eigenvalues in an increasing order of {£ £ } in Q. We will show in Section 7 that 

— /ig| < CkS if f2 is C 1,1 (or convex in the case m = 1), and — yUpl < Cfc )0 .e| ln(e)|2+°" 
for any cr > if Q is Lipschitz. The same holds for Dirichlet and Steklov eigenvalues. To 
the best of the authors' knowledge, only results for Dirichlet eigenvalues in smooth domains 
[T0| [15] and Neumann eigenvalues in a two-dimensional curvilinear convex polygon [16] were 
previously known (see [TTJ, [18j [191 120] for related homogenized eigenvalue problems). 

Finally, in Section 8, we prove several weighted L 2 potential estimates, which are used 
in earlier sections, for the operators {C £ }. Our proofs use asymptotic estimates of the 
fundamental solutions for C £ in [3] as well as some classical results from harmonic analysis. 

The summation convention is used throughout this paper. Unless otherwise stated, we 
always assume that A e A(/x, A,r), A* = A, and Q is a bounded Lipschitz domain in R d , 
d > 2. Without loss of generality we will also assume that diam(fi) = 1. We will use C 
and c to denote positive constants that depend at most on d, m, /i, A, r and the Lipschitz 
character of Q. 



2 Uniform regularity estimates 

In this section we recall several uniform regularity estimates for {£ e }, on which the proofs 
of our main results rely. We also give definitions of the non-tangential maximal function and 
radial maximal operator Ai. 

Let u e be a weak solution of C e {u e ) = in Q. Then if B(x, 2r) C Q, 

\Vu £ (x)\<-^ [ \u E (y)\dy. (2.1) 

r JB(x,r) 

This uniform gradient estimate was proved in [1] (the symmetry condition A* = A is not 
needed for this). Let T £ (x, y) = (r° /3 (x, y)) denote the fundamental solution matrix for C £ in 
]R d , with pole at y. It follows from the gradient estimate (12.11) that \T £ (x, y)\ < C\x — y\ 2 ~ d , 
\V x T e (x,y)\ + \V y T e (x,y)\ < C\x ~ y\^ d and | V x V y T e (x, y)\ < C\x - y\~ d (see [3]). 

For a function u in a bounded Lipschitz domain Q, the non-tangential maximal function 
(u)* on dQ is defined by 

(u)*(Q) =sup{\u{x)\ : x G D, and \x - Q\ < C dist{x,dQ)} , (2.2) 
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where Co, depending on d and the Lipschitz character of f2, is sufficiently large. 

Theorem 2.1. Let f G L 2 (dQ) and u £ be the unique solution of the L 2 Dirichlet problem: 
C £ {u £ ) = infl, u £ = f non-tangentially on dQ and (u £ )* G L 2 (dQ). Then 

II (ue)*\\L*(dn) + 11^11^1/3(0) + jy \Vu £ (x)\ 2 5(x)dx^j < C\\f\\ L 2 {dnh (2.3) 

where d~(x) = dist(x, <9fi). Furthermore, if f G ^(dQ), i/ie solution satisfies the estimate 
\\(Vu E )*\\ L 2 m < C\\f\\ H i m . 

Proof. The non-tangential maximal function estimate ||(« e )*||L 2 (9f2) — ||/||l 2 (9Q) i n Lipschitz 
domains was proved in [5] for m = 1 and in [TS] for m > 1. In the case of smooth domains, 
the estimate was obtained earlier in [21 1]- The proof of || (Vu £ )*\\ L 2^ a ^ < C||/||^i(an) may 
be found in [12] for m = 1 and in [13] for m > 1. 

It was also proved in [13] that the solution of the Dirichlet problem with boundary 
data / in L 2 (dQ) is given by a double layer potential V £ (g £ ), where the density g £ satisfies 
I belli 2 (an) < C||/||i 2 (9n)- This, together with Proposition 18.5} gives the square function 
estimate in (12. 3p . 



[J \Vu £ (x)\ 2 8{x)dx\ <C\\g £ \\ L 2 {dn) <C\\f\\ L 2 {m) . 



Finally, the estimate ||w£||#i/2(n) — C\\f\\L2(dn) follows from the square function estimate by 
real interpolation (see e.g. [9] pp. 181-182]). □ 

The next theorem was proved in [13] (the case m = 1 was obtained in [12]). We refer the 
reader to [HI [TJl d3] for references on L p boundary value problems in Lipschitz domains in 
non-homogenized settings. 

Theorem 2.2. Let g G L 2 (d£l) with J m g = 0. Let u £ G H\Q) be the unique (up to an 
additive constant) weak solution of the L 2 Neumann problem: C £ (u £ ) = in ; ^ = g on 



tan 

AT^oirn nrtrt mrnhlprn- i (it \ — PI on O 

due 

dtt. Then ||(Vu e )*IU 2 (an) < C||3>|| L 2 (an) . 



The radial maximal operator. Given a bounded Lipschitz domain Q, one may construct a 
continuous family {Q t , —c < t < c} of Lipschitz domains with uniform Lipschitz characters 
such that Qq — ^ an d fit C f2 s for t < s. We may further assume that there exist homeo- 
morphisms A t : dtt -» 0Cl t such that A (Q) = Q, \A t (Q) - A S (P)\ ~ \t - s\ + \P - Q\ and 
|A S (Q) — A t (Q)\ < C odist (A s (Q) , d fl t ) for any t < s (see e.g. [22]). For a function u in Q, 
the radial maximal function A4(u) on dfl is defined by 

M{u){Q) = sup{|u(A t (Q))| : -c<t<0}. (2.4) 

Observe that Ai(u)(Q) < (u)*(Q) and if S C f2 is a surface near dfl and obtained from 
dfl by a bi-Lipschitz map, then ||w|| L 2( 5 ) < C 1 1 (n) 1 1 (qq) . Also, note that ||w||l 2 (q) + 
\\M{u)\\ L 2 (dn) < C||(w)*||i 2 (9n)) an d the converse holds if u satisfies the interior L°° estimate 
||w|U°°(B) < C|2_B| _1 ||u||ii( 2j B) for any 2B C Q. In particular, if C £ {u £ ) = in Q, then 

|| (U e )* || L 2 (3fi) ~ |K|U 2 (fi) + \\M{u £ )\\ L 2 {m) . 
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3 Homogenization of elliptic systems 

Let C £ = — div(A(x / e)W) with A(y) satisfying (jl.2p ~ fll.3p . The matrix of correctors xiv) 

a/3 



(Xj(y)) f° r {£-e} is defined by the following cell problem: 



dy 

,0 

3 



a/3 , aj 
a ij + a ik 



d 

dy k 



in 



a 



Xj^iy) is periodic with respect to Z d , 
J xfdy = 0, 



(3.1) 



for each 1 < j ' < d and 1 < /3 < m, where Y = [0, l) c 
is given by Cq = — div(AV), where A 

~a/3 _ 



(a*f) and 



1 /Z d . The homogenized operator 



aB i cry 



dyk 



xf 



dy 



(3.2) 



(see i). 



Lemma 3.1. Lei F = (F 1 ,...,F d ) G L 2 (Y). Suppose that J y Fj dy = and div(F) = 0. 
Then there exist Wij G H l (Y) such that Wij 



-Wji and Fj 



duii. 
dyi 



Proof. Let fj G H 2 (Y) be the solution to the cell problem: Afj = Fj in Y, fj is periodic 
with respect to Z d and j Y fj dy = 0. Since div(F) = 0, we may deduce that is constant. 
From this it is easy to see that 

df, df\ 



dyi dyj 



has the desired properties. 
Let 







<(y) = ^-ojf(y)--C(y)^-( X f). 



□ 



(3.3) 



It follows from (jH2]) and (ED) that 

J$<*fdy = and A($°f)=0 in R d . 

Hence we may apply Lemma 1331 to (y) (with a,(3,j fixed). This gives G Lf 1 (F), 
where 1 < i, j, k < d and 1 < a, < m, with the property that 



a/3 
ij 



d 

dy k 



«g} ^d *g 



(3.4) 



Furthermore, it follows from the proof of Lemma 13.11 that if x £ W /1,P (V) for some p > d, 
then ^ G L°°(F). 

The next lemma is more or less known (see e.g. [TOl Chapter 1] for the case m = 1 and 
v e — u o)- We provide the proof for the sake of completeness. 
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Lemma 3.2. Let u a e G H 1 ^), G H 2 {9) and 



" ( x / £ )^ 

where 1 < a < m. Suppose that C E {u £ ) = Cq(v e ) in Q. Then 



w~(x) = <(*) - vf(x) - ex a k p (x/e)-^, (3.5) 



(3.6) 



where 



K; k (y) = ^ k (y) + ( hf(y)xr(y) (3.7) 

and ^Ji k (y) is given in \3.4\ )- 

Proof. It follows from the assumption C £ {u £ ) = Co(v £ ) that 

^k{< {x/ ^ {xM k^k} 
- ~k + i ' 

where the periodic function $f fc 7 (y) is given by (13. 3p . Using the first equation in (13. 4p . we 
obtain 

(^* = ^{w<} 2 (38) 

+ £ |-{(*S (l/£) + < ( ^ 7(l/£) )^}- 

By the second equation in (13. 4p . the first term in the right hand side of (I3.8P is zero. This 
gives the equation (13.61) . □ 

Remark 3.3. Under the assumption A G A(/i, A,r), it is known that Vx is Holder con- 

'°^ k is Holder continuous. In particular, ^f k , 

T. 

Fix F G L 2 (fi) and / G H l ' 2 {dty. Let w £ , m G H\Q) solve 

£ £ ( Me )=F inO, ^ (Co(u ) = F in S7, 



tinuous. This implies that V^°f k is Holder continuous. In particular, ^f k , b^ k G L°°(Y). 
Furthermore, ||^"^||oo + ll^^lloo is bounded by a constant depending only on m, d, /i, A and 



u £ = f on d£l, ' 1 Uq = / on <9fi, 

respectively. 

Theorem 3.4. Lei be a bounded Lipschitz domain. Suppose that A G A(/i, A, r) and 
A* = A. Assume further that uq G H 2 (Q). Then 

\\M(w £ )\\ L 2 {m) + \\w £ \\ H i/2 {Q) + |y |Vw £ (x)| 2 5(x) da; | < C£||w ||jf2 (n) , (3.10) 
where w £ (x) = u £ (x) — Mq (x) — ex^{x/e)^- and 5(x) = dist(x, dfl). 



Observe that 



\\U £ - U \\ L 2 {n) + \\M(u e - U )\\ L 2 (dn) 

< \\w e \\&({i) + \\M(w £ )\\ L 2 {dn) + Ce {\\Vu \\ L2{n) + ||A^(Vn )|| L2(en) } (3.11) 

< ||^e|U 2 (n) + \\M(w £ )\\ L 2 {m) + Ce\\u \\ H 2 {n) , 

where we have used the fact that ||A'l(VMo)||L 2 (an) < C|l M o||_H" 2 (n)) which follows from the 
estimate (18.151) . 

As a corollary of Theorem I3.4[ we obtain the following. 



Corollary 3.5. Under the same assumptions as in Theorem 3.4, we have 
\W - u o\\^(n) + \\M{u £ - u )\\ L 2 {dn) < Ce\\u \\ H 2 {n) . 



Proof of Theorem\3.4\ We first observe that by (13.61) . w £ satisfies 



{C £ (w £ )} a = e- 9 



d 2 ul 



dxi \ ljk dxjdx k 



in Q, 
on d£l. 



(3.12) 



(3.13) 



Let w = 9 £ + z £ , where 



{C £ (6 £ )} a = in ft, 



du 



dxk 



(3.14) 



and 



To estimate 6 £ , we apply Theorem 12.11 to obtain 



M(6 £ )\\ L 2 (m) + ||0 e || H i /2(n) + i / \V6 £ (x)\ 2 5(x) dx \ < Ce\\Vu \\ L 2 



1/2 



(3.15) 



(3.16) 



Since tf*? k G L°°(Y), by the usual energy estimates, we have H^H^i^) < Ce\\ V 2 wo||,L 2 (n)- 
Thus, 



\\M(z £ )\\ L 2 {m) + \\z £ \\ H i/2 (n) + |y \Vz £ (x)\ 2 5(x)dx^ < Ce\\V 2 u \\ L 2 {n) . (3.17) 
Since || VMo||L 2 (ar2) < C||wo||,tf 2 (f2), the desired estimate (13. lOj) follows from (I3.16P and (I3.17p . 



This completes the proof of Theorem 13.41 



□ 
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Remark 3.6. Let Q be a C 1,a domain in M. d . As we mentioned in the Introduction, the 
estimate \\u £ — Wo||L 2 (n) < C£|l w o||i/ 2 (n) was proved in [T5], using the estimates for the L 2 
Dirichlet problem in [2J [Tj . Let 9 £ and z £ be given by ( 13.1 4p and (I3.15P respectively. It follows 
from pQ Theorem 3] that H^H^oo^) < Ce\\ Vuo\\L°°(dn)- In view of (13 . 1 5[) we have 



\z £ (x)\ <Ce [ \V y G £ (x,y)\ \V 2 u (y)\dy, (3.18) 



where G £ (x,y) denotes the Greeen function for C £ in Q. By [T] we have \V y G E (x, y)\ < 
C\x — y\ 1 ~ d . It follows from (13.181) and Holder inequality that ||2 £ ||.L°°(n) < C p e\\ V 2 uo\\lp(u) 
for any p > d. This gives 

\\ue - «o|U°°(n) < C p e\\u \\ W 2, P{n) foranyp>d, (3.19) 

where we also used the Sobolev imbedding || VwoII^q) < Cpll M o|| w 2 'P(n) for P > d. 

4 Dirichlet boundary condition 

Let Q be a bounded Lipschitz domain. Let Q £ D Q and A e : dQ — > dQ £ be defined as in 
Section 2. Given / G H l (d£l) and F G L 2 (fi), let u £ G -H 1 ^) be the weak solution of 

£o(^ £ ) = ^ in fi e> C4 x) 

w e = f £ on <9fi £ , 

where F = F in f2 and zero otherwise, and f £ (Q) = /(A" 1 ^)) for Q G <9f2 £ . The goal of 
this section is to prove the following. 

Theorem 4.1. Let Q be a bounded Lipschitz domain. Suppose that A G A(/i, A,r) and 
A* = A. Let 



w 



»(*) = <(*) - <(*) - e X f (x/e)p-, (4.2) 



where v £ is given by ( |^.i| ). T/ien, if < e < (1/2), 

| L 2 (n) <Ce|ln(e)| a {||F|| L 2 (n) + 11/11^1^)}, for any a > 1/2, (4.3) 



||.M(OIU W < Cerl \n(e)\ a {\\F\\ L 2 (U) + ||/|| H i(8n)} , for any a > 3/2, (4.4) 

and 

\\we\\w/*(n) + |y \Vw £ (x)\ 2 5(x)dxj < Ce\\ne\ {||F|| Z 2 (Q) + ||/||tfi ( an)} , (4.5) 

where C depends only on \l, \, r, d, m, a and Q. 

As a corollary we obtain the following convergence rates of u £ to uq in L 2 . 
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Corollary 4.2. Under the same conditions as in Theorem 4-1, we have 

\\u £ -u \\ L 2 (n) < CE\\n(s)\^{\\F\\ L2{n) + \\f\\ Hlm) }, (4.6) 

\\M(u £ - u )\\ L , m < Ce\ \ne\* + ° {\\F\\ LHn) + ||/|| H 1 {dQ)} , (4.7) 

for any a > 0. 

Without loss of generality we shall assume that H-FUl 2 ^) + 11/ 1 1 if 1 (an) = 1 in the rest of 
this section. We begin with an estimate on Vt> £ . 

Lemma 4.3. Let v £ be defined by ^4-l\) . Then 

\\M £ {Vv £ )\\ L 2 {mE) + \\Vv e \\ m ,2 m + |jf \V 2 v £ (x)\ 2 5 £ (x)dxj < C, 

where 5 £ (x) = dist(x,dVt £ ) and M £ (v £ )(A £ (Q)) = sup {\v £ (A s (Q))\ : — c < t < e}. 

Proof. Let G = T Q * F in R d , where T (x) is the matrix of fundamental solutions for the 
operator £ , with pole at the origin. Clearly, ||G||h 2 (k <1 ) — C||-^1U 2 (n)- This implies that 
\\M e (VG)\\ LHdQe) + \\G\\ H i (eQe) < C\\F\\ L 2 (n) . 

Next, we note that C (v £ — G) = in Q £ and v £ — G = f £ — G on dfl £ . Hence, by Theorem 
12.11 (see [7] for operators with constant coefficients), 



||(V(^ - G)Y\\ LHdns) + \\V(v £ - G)\\ m/Hns) + | jf |V 2 '- ^ 



(v £ -G)\%(x)dx 



< C\\f £ - G\\ H i( d n £ ) < C. 

It follows that 

1/2 



M £ {Vv £ )\\ L 2 {dns) + \\Vv e \\ H i /a({le) + <^ / \V 2 v £ {x)\ 2 5 £ {x)dx 

<C+\\M £ (VG)\\ L 2 me) + \\VG\\ H , mUe ) + \V 2 G(x)\ 2 5 £ (x)dx} 



1/2 



< c. 



□ 



Remark 4.4. By Lemma [4.31 we have ||Vu e ||i2m) + ||A^(Vu e )||i,2(an) < C. It follows that 

IK - Uo||x2(n) < + Ce, 

\\M(u s - u )\\ L 2 m) < \\M(w £ )\\ L 2 (m) +Ce. 



This, together with Theorem I4.1[ gives the estimates in Corollary [4] 

Lemma 4.5. Let f £ (Q) = /(A~ 1 (Q)) and v £ be defined by Then \\f — u £ ||L 2 (an) < Ce 

and ^ ^ 

||(v e -ito)*||£a(8n) + ||«e-«o||Hi/a(n) + |y |V(u e -Uo)| 2 £(ac)da:J < Ce. 
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Proof. Note that for Q G 

\f(Q) - v £ (Q)\ = \v £ (Ae(Q)) - v £ (Q)\ < CeM £ {Vv £ ){K{Q))- 

This gives ||/ — w £ ||l 2 (<9q) < CE\\M. £ {yv £ )\\i,2(dn. e ) < Ce, where the last inequality follows 
from Lemma 14.31 Since Cq{v £ — u ) = in Q and v £ — u = v £ — f on dQ, we may apply 
Theorem 12.11 (for the case of constant coefficients) to obtain 

■ 1/2 

(v e - u y\\ L 2 {m) + \\v £ - «o 11^1/3(0) + <j I |V(v e - u )\ 2 5(x) dx 

<C\\v £ -f\\ L 2 m <Ce. 
This completes the proof. □ 
Let <f) a (t) = {ln(i + e a )} a . 

Lemma 4.6. Let W £ G i? 1 (f2) 6e a solution of C £ (W £ ) = div(h) in f2 and PV £ = g on d£l 
for some h G L 2 (fi) and g G H 1 ^^). Tnen 

ITOU 2 ^) < C||^|| x2(9 n) +C *\J \h{x)\ 2 5{x)(j) a {5{x)) dxl for any a > 1, (4.9) 

/■ 1/2 (- /■ 1/2 

/ |V^ £ (x)| 2 5(a;) da; < C\\g\\ L * m +C{ \h{x)\H{x)<p 2 {5{x)) dx , (4.10) 
and 

1/2 



\\M(W £ )\\ L 2 {dn) < C\\g\\ L 2 m + C a |^ |/i(x)| 2 5(x)0 a (5(x))dxj (4.11) 

/or any a > 3. 

Proof. Let ft, = (/if) and 

#f (*) = - jT ^ {Tf (x, y)} fcf (y) dy, (4.12) 

where r e (x,?/) = (r" /3 (x,y)) is the matrix of fundamental solutions for C £ in R d , with pole 
at y. Note that C £ (W £ — H £ ) = in Q and W e — H £ = g — H £ on <9fi. It follows by Theorem 
EH that 

f r 1V2 
||W £ - # £ |b(n) + ||(W e - J ff e )*|| za(ai) + ^ / |V(iy e -if £ )| 2 5(a;)da; 



< C||<?IU 2 (an) + C\\H £ \\ L 2( dn y 

Hence, 

HW^IU^n) < C { ||o||L 2 (9n) + ll-ffelU 2 ^)} + ||-He|U 2 (n), 



\\M(W £ )\\ L 2 im) < C{\\g\\ LHdn) + \\M(H £ )\\ LHdn) } , 
/ \VW £ \ 2 S(x)dx) <C{\\g\\ LHdn) + \\H £ \\ LHdn) } + { [ \VH £ \ 2 S(x)dx 
The desired estimates now follow from Propositions 18.11 18.21 18.31 and 18.41 □ 
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We are in a position to give the proof of Theorem 14. 1[ 
Proof of Theorem J^.l, Let 

f) v P 

fi? = <(*) - t£(x) - exfix/e)-^ = < - (« e - u ) (4.13) 
in f2. In view of Lemma 14.51 it suffices to show that w e satisfies the estimates in Theorem 

El 

To this end we first observe that by Lemma | 



Let w £ = 9 £ + z £ , where C £ {9 £ ) = in fi, 9 £ = w £ on dCl, and z £ satisfies 

To estimate 9 £ , we apply Theorem 12.11 to obtain 

ll^l|L 2 (n) + ||^(^)IU 2 (^) + {^|V^| 2 5(a:)dx| ' <C\\9 e \\ L , m) 

= C\\w £ \\ L 2 {m) < C{\\f - v £ \\ L2{dn) + e\\Vv £ \\ L 2 {m) } < Ce, 

where the last inequality follows from Lemmas 14.31 and 14.51 

Finally, we use Lemma [4.61 to handle z £ . In particular, this gives 

WzeWvp) < Ce I jf \V 2 v £ \ 2 5(x)(/> a (5(x)) da: J ' . (4.14) 

for any a > 1. Note that a (£) is decreasing and t(p a (t) is increasing on (0, oo) for any a > 0. 
Hence, for any i6l] and < e < cq, 

6(x)(f) a (x) < 5 £ {x)4> a {5 £ {x)) < 5 £ {x)4>a{e/C) < CS £ (x)\ \n(e)\ a , 

where 5 £ (x) = dist(x, dft £ ). In view of (I4.14p we obtain 

\V 2 v £ \ 2 5{x)(t)a{5{x)) dx < Ce\ \n(e)\ a [ \V 2 v £ \ 2 5 £ {x) dx 

Jn 

< Ce\\n(e)\ a [ \V 2 v £ \ 2 5, 

< Ce\ln{e)\ a , 



(4.15) 



for any a > 1, where the last inequality follows from Lemma 14.31 Thus ||z £ ||l 2 (q) < 
Ce\ ln(e)| a//2 for ant a > 1. This, together with the estimates of 9 £ and v £ — u in L 2 (Q), 
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gives (14.31) . Estimates ( 14. 4p and ( 14. 5 p follow from Lemma [4.61 in the same manner. We omit 
the details. □ 

Proof of Theorem li.il Estimate (11.51) is given in Corollary 13.51 and estimate (jl.6p in 
Corollary OJ □ 



5 Neumann boundary condition, part I 

Fix F G L 2 (Q) and g G L 2 (dfl). Suppose that f Q F + j m g = 0. Let u £ , u G H 1 ^) solve 

( C £ (u £ ) = F intt, ( C (u ) = F inn, 

\fe = g on dn, and \$* = g on dn, V- 1 * 

respectively. Recall that 



where n = • ■ ■ , rid) denotes the outward unit normal to dfl. 

Lemma 5.1. Letw? = v?-vf- exf(x/e)^-, where u £ G H 1 ^) and v £ G H 2 (n). Then 

1 - n^xja^r (x/e)- n^xja^ 



dv £ ) »w ,J v / > q x . .V / *J ^ 

l{ n * (a;) ^ 



d 2 v 1 

em(x)bJ k (x/e)- 



dxjdxk 

where ^^{y) and b^ k {y) are the same as in Lemma WM. 
Proof. A direct computation shows that 



afif I \ dw e I \ ap, , \ du e / \^ap^e 



+ TH{x)*%{x/e)?j£- (5.4) 
- en i (x)a^(x/e)xt'{x/e)- 



dxjdxk 

where ^ k {y) is defined by (13.31) . By (13.41) . we obtain 



<9x fc 9a;j [ J dx k 



d 2 v"< 

-eni(x)^ k (x/e)- 



dxjdxk 
d 2 v~ 1 

-en^Ux/e)- ; 



dxjdxk 
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(5.5) 



Equation (15. 3p now follows from ( 15. 4 p and (15. 5p . □ 
Theorem 5.2. Let Q be a bounded Lipschitz domain. Suppose that A G A(/i, A, r) and 



A* = A. Let (u £ , Wo) ^ e a solution of ( 15.jp u>i£/i 
m G # 2 (ft). T/ien 



90 M e 



Idn u ° = ®- Assume further that 



M(w e )\\ L 2 {dn) + \\w £ \\ H i/2 (n) + < / |Vw e (x)| [• <Ce\\u \\ H 2 



1/2 



2(C), 



(5.6) 



where = u £ (x) — Uq(x) — e~x^f {x / e)^- . 

As in the case of Dirichlet boundary conditions, Theorem 15.21 gives the following con- 
vergence rate of u £ to uq in L 2 . As we mentioned in the Introduction, the estimate ||it E — 
w o|| j l 2 (S7) — Ce\\uo\\h 2 (q) was proved in [16] when ft is a curvilinear convex domain in R 2 . 

Corollary 5.3. Under the same assumptions as in Theorem \5.!Q we have 

IK - u \\ L 2 {n) + \\M{u £ - u )||i»(an) < Ce\\u \\ H 2 {n) . (5.7) 

Proof of Theorem I5.H 

In view of Lemmas 13.21 and 15.11 we may write w e = 9 £ + z £ + p, where 

( C £ (6 £ ) = in a 



(Me 



- \ nAx 



_d_ 

dx„ 



on <9ft, 



(5i 



9 £ G F x (ft) and [ 9 £ = 0, 
Jan 



and 



a 2 < 1 

|^ fcV ' ~' dxjdx k J 
-eni(a;)6g(a;/e)- 

z £ G tf^ft) and / 

Jn 

1 



dz £ 
dv f 



d 2 ul 

dxjdxk 
0. 



in ft, 
on c?ri 



(5.9) 



P 



\dtt\ 



(w e - z £ 



is a constant. It follows from the energy estimates that ||^ e ||jfi(n) < C^llwoll^^). Also note 
that 



p\<C \z £ \ +Ce \Vu \ < Ce\\uo\\ H 2 



mi 



»(n), 



where we have used the condition f. 



on Ue 



Idn u ° = ^" Thus it remains only to estimate 9 £ . 



To this end we use a duality argument and consider the L 2 Neumann problem 

( C £ (Q £ ) = in ft, 

= /i on 9ft, 



96, 



<9z/. 



(5.10) 



6 £ G if^ft) and / 6 e = 0, 
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where h G L 2 (dQ) and J an h = 0. It follows from integration by parts that 



Oe-h 



Ml 



,,<> dis £ * 



d d 
Tii— n« 



^7 



(5.11) 



< Ce||Ve e || L 2 (a n)||Vuol|L2 



where we have used the fact that tangential derivative for 1 < i, j < d. In 

view of Theorem 12.21 we have || VO e ||i2( 9n ) < C||/i||£2( 9 n). Hence, by ( 15. lip and duality, we 
obtain ||0 e ||.L2(an) < C£||Vtto|U 2 (an)- Here we also use the fact f dn 9 £ = 0. 

Finally, we use the estimates for the L 2 Dirichlet problem in Theorem 12.11 to see that 

Peh^U) + \\M(9 e )\\v(n) + \V6 £ (X)\ 2 5(X) dx} ' 

< C||^I|l 2 (9c) < Ce\\\7u \\ L 2( dn j < Ce||«o||.ff2( n ). 
This, together with the estimates of z £ and p, completes the proof of Theorem 15.21 □ 



Remark 5.4. The estimates in Theorem 15.21 and Corollary 15.31 also hold under the condition 
In Ue = In u ° = ®- ^ n ^ n * s case ^ ne cons t an t P is given by p = ^ Jq(w £ — 9 £ ), and we have 

\p\ <Ce |V«o| + C \6 £ \ < Ce\\u \\ H 2 {n) . 
Jn Jn 

This will be used in the proof of the error estimate for the Neumann eigenvalues for L £ . 



6 Neumann boundary condition, part II 

In this section we extend the results on convergence rates in Section 4 to the case of Neumann 
boundary conditions. 

Construction of the first-order term. Fix F G L 2 (Q) and g G L 2 (dQ) such that J Q F + 
f dn g = 0. Let (u £ ,Uo) be the solution of (15. ip with f m u £ = J dn uo = 0. Consider 

< = <-<- £X f(z/e)^, (6.1) 

where v £ is the solution of (14. ip in Q £ with / given by u \q^ and f £ (Q) = /(A~ 1 (Q)). Note 
that by Theorem 12.21 (for operators with constant coefficients), 

11/11^1(00) < C {\\g\\ L 2 (m) + ||F||t2 (9n) } . (6.2) 

Theorem 6.1. Let Q be a bounded Lipschitz domain. Suppose that A G A(/x, A,r) and 
A* = A. Let w £ be defined by (GOP- Then, if < e < (1/2), 

\\w £ \\ L \dn) + IKIU^n) < Ce\ \n(e)\ a {\\F\\ L 2 {n) + ||^||£2 (an) } , for any a > 1/2, (6.3) 
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\\M{w £ )\\ L 2 {m) < Ce\ \n(e)\ a {\\F\\ L 2 {n) + \\g\\ L 3 (gn) } , for any a > 3/2, (6.4) 



and 



\w F 



HVa(n) + |y |Vw £ (x)| 2 5(x)rfx| < Ge|lne| {||F||jy«(n) + ||<7||j^(an)}, 



(6.5) 



where C depends only on \l, \, r, d, m, a and Q. 
Observe that 

||Vv e ||£a ( n) + \\M(Vv £ )\\ LHm) < C {\\F\\ L 2 {n) + ||/|| H i (S n)} 

< C {ll-^IU 2 (Q) + lblU 2 (sn)} • 



Thus, as a corollary of Theorem 16.14 we obtain the following convergence rates of u £ to uo 
mL 2 . 

Corollary 6.2. Under the same conditions as in Theorem \6.1l we have 

\\ue - Holly (an) + IK - «o 11^(0) < Ce\ ln(e)|^ +CT { \\F\\ L 2 {n) + H^Hx^can) } , (6-6) 
I |.M (u 6 - u )\\ L 2 {m) < Ce\ lne|i +a {||F|| L2(n) + ||#|| L 2 (8 n)} , 

for any o > 0. 



(6.7) 



Without loss of generality we will assume that H-FH^/q) + || i»ra«) < 1 in the rest of 
this section. We remark that because of (16. 2p . the estimates of Vv £ in Lemmas 14.31 and 14.51 
continue to hold. 

Proof of Theorem \6.1\ We proceed as in the case of Dirichlet condition and write 



< = < + {< -Mq} and w« = u«(x) -v?(x) -exf(x/e) 



dx k 



The desired estimates for v £ — Uq follow directly from Lemmas I4.3ti4.5l and (16. 2p . 

Next we let w £ = 9 £ + z £ + p, where 



e\ w e ) 







in Q, 



r ft 



_d_ 



(6.8) 



G H\Q) and / 9 £ = 



{CsCzs)T = e^\b^ k {x/e) 
dz £ \ a 



d 2 v] 
dxjdxk 

d 2 v] 
dxjdxk 



in Q, 



on dfl 



(6.9) 



z £ E H\tt) and / : 

JdQ 



z s = 0, 
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and 



is a constant. Note that 



\dn\ 



an 



\dn\ J m \dQ\ J dn 



\p\ < Ce\\Vv £ \\ L 2 (m) + C\\v £ - u \\ L 2 {m) < Ce, 

where we have used the fact J dn u £ = J an u = as well as Lemmas 14.31 and 14.51 
By a duality argument similar to that in the proof of Theorem 15.21 we have 

||#£||L 2 (0n) — Ce\\ Vf e ||L 2 (on) < Ce. 
It then follows from the estimates for the L 2 Dirichlet problem in Theorem 12.11 that 

1/2 

\\M(e £ )\\ L 2 idn) +\\e £ \\ H1/2(n) + Ujve £ \ 2 5(x) dx ' 

< C||^elU 2 (an) < Ce. 

The estimates of z £ also relies on a duality estimate. Indeed, let O e G be the 

solution of (15.1 Op with h G L 2 (dQ) and f Qn h = 0. It follows from integration by parts that 



z £ ■ h 



mi 



oe £ 



an du £ 



(6.10) 



d 2 v~< dO a 

&/e) 



dxjdxk dxi 



By the Cauchy inequality this gives 



z £ ■ h 



1/2 

2„ |2.c, ...... ( J / l V Q |2 



<e{ j \V 2 v £ \ 2 5(x)<p a (5(x))dx 



dx 



5{x)<j> a {5{x)) 



1/2 



(6-H) 



Observe that if a > 1, 



1/2 



5(x)0 a (5(x)) 
Hence, by (16.1 ip and duality, we obtain 



z E \\L^au) <Ce< / |V v £ \ 8(x)(j) a (S(x)) dx 



1/2 



(6.12) 



for any a > 1. With (I6.12p at our disposal we may apply Lemma H~6l to obtain 

|*e||2/»(n) ' " 



<Ce|y |VV| 2 W0a(<Kx))dxj < Ce\ ln(e)| a/2 , 
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where we have used (I4.15P for the last inequality This, together with estimates of 9 £ , p and 
v £ — uo, gives ( 16.31) . Estimates ( 16. 4 p and ( 16. 5 p follow from Lemma 14. 6 1 and (16.121) in the same 
manner. This completes the proof of Theorem 16. 1[ □ 

Remark 6.3. The estimates in Theorem 16.11 and Corollary 16.21 continue to hold under the 
condition f n u £ = f Q Uo = 0. In this case one has 

Hence, 

|p| < Ce||Vu e ||L3(n) + C\\v £ - u \\ L ^n) + C\\9 £ \\ L 2 {n) + C\\z £ \\ L 2 {n) 
< Ce\\n{e)\ a/2 . 
for any a > 1. The rest of the proof is the same. 

Proof of Theorem [7TE Estimate (II. 5p for the Neumann boundary conditions is given in 
Corollary 15.31 and estimate (II. 7\i in Corollary 16.21 □ 



7 Convergence rates for eigenvalues 

In this section we study the convergence rates for Dirichlet, Neumann, and Steklov eigenval- 
ues associated with {C £ }. Our approach relies on the following theorem, whose proof may 
be found in [TOj pp.338-345]. 

Theorem 7.1. Let {T £ ,e > 0} be a family of bounded, positive, self-adjoint, compact op- 
erators on a Hilbert space 7i. Suppose that (1) \\T £ \\ < C and \\T e f — T f\\ — > as e — > 
for all f G H; (2) {T £ f £ ,e > 0} is pre-compact in H, whenever {f £ } is bounded in 7i. Let 
{u £ } be an orthornomal basis of Ti consisting of eigenvectors of T £ with the corresponding 
eigenvalues {fi £ } in a decreasing order, 

/4>^ 2 >--->^>--->o. 

Then p £ > Ck > 0, and if e > is sufficiently small, 

\l4 ~ Pol ^ 2su P {\\ T eU ~ T u\\ : u e N(/Iq,T ) and \\u\\ = l}, 
where N(fiQ,To) is the eigenspace of To associated with eigenvalue Hq. 

Dirichlet eigenvalues. Given / £ H = L 2 (Q), let Tf(f) = u £ G Hl(VL) be the weak solution 
of C £ (u £ ) = f in Q. It is easy to see that {T^,e > 0} satisfies the assumptions in Theorem 
17.11 Recall that X £ is a Dirichlet eigenvalue for C £ in Q if there exists a nonzero u £ G (Q) 
such that C £ {u £ ) = \ £ u £ in Q. Let {A^} denote the sequence of Dirichlet eigenvalues in an 
increasing order for C £ . Then {(A^) -1 } is the sequence of eigenvalues in a decreasing order 
for T® on L 2 (Q). It follows from Theorem 17. II that if e is sufficiently small, 

< 2 sup {\\u £ - tto||i2 ( n)}, (7.1) 



1 1 

aF" Ag 
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where the supremum is taken over all u E , u G H^ft) with the property that Aq ||ii ||i 2 (r2) = 1 
and C £ {u £ ) = £o(u ) = XqU in Q. Note that if Q is C 1 ' 1 (or convex in the case m = 1), 
then 1 1 uq \\h 2 (v.) < Cfc by the standard regularity theory for second-order elliptic systems with 
constant coefficients. Also, by Theorem 17.11 we see that X k < Ck- Hence, we may deduce 
from (17.11) and Corollary 13.51 that |A* — A§ | < c^e (see e.g. [TQl p. 347]). However, if is a 
general Lipschitz domain, uq G H 2 (Q) no longer holds. Nevertheless, Corollary 14.21 gives us 
the following. 

Theorem 7.2. Let < e < (1/2) and {X £ } be the sequence of Dirichlet eigenvalues in an 
increasing order of C £ in a bounded Lipschitz domain. Then for any o > 0, 

|A*- Aj | < ce\ln(e)\? +a , 

where c depends on k and a, but not e. 

Neumann eigenvalues. Given / G L 2 (Q) with J ' / = 0, let T^(f) = u £ G H l (Q) be the 
weak solution of the Neumann problem: C £ (u £ ) = f in fi, |^ = on <9f2 and J n u £ = 0. 
Again, it is easy to verify that the family of operators {T^} on the Hilbert space {/ G 
L 2 (Q) : J Q f — 0} satisfies the assumptions of Theorem 17.11 Recall that p £ is a Neumann 
eigenvalue for C £ in Q if there exists a nonzero u £ G if x (f2) such that C £ (u £ ) = p £ u £ in f2 and 
|^ = on dVl. Let {p^} be the sequence of nonzero Neumann eigenvalues in an increasing 
order for C £ in Q. Then {(p^)^ 1 } is the sequence of eigenvalues in a decreasing order for . 
Thus, in view of Theorems 17.11 and Remarks 15.41 and 16.31 we obtain the following. 

Theorem 7.3. Let < e < (1/2) and {p £ } denote the sequence of nonzero Neumann 
eigenvalues in an increasing order for C £ in a bounded Lipschitz domain Q. Then for any 
a>0, 

\p k £ -p k \<ce\Me)\^, 

where c depends on k and a, but not e. Furthermore, the estimate \p £ — p^\ < CkS holds if 
Q is C 1 ' 1 (or convex in the case m = 1). 

Steklov eigenvalues. We say s £ is a Steklov eigenvalue for C £ in Q if there exists a nonzero 
u £ G H 1 ^) such that C £ {u £ ) = in and §^ = s^dQl' 1 ^ on dQ. Note that s^^l" 1 is 
also an eigenvalue of the Dirichlet-to-Neumann map associated with C £ . Given g G L 2 (dQ) 
with f Qn g = 0, let S E (g) = u £ \dn, where u £ is the weak solution to the L 2 Neumann problem: 
C £ (u £ ) = in f2, |^ = g on dil and J dn u £ = 0. It is not hard to verify that the family 
of operators {S £ } on the Hilbert space {g G L 2 (dQ) : J dn g = 0} satisfies the assumptions 
in Theorem 17.11 Consequently, the L 2 (dQ) convergence estimates in Corollaries 15.31 and 16.21 
give the following. 

Theorem 7.4. Let < s < (1/2) and {s £ } denote the sequence of nonzero Steklov eigen- 
values in an increasing order for C £ in a bounded Lipschitz domain Q. Then for any a > 0, 

\s k e -4\ < ce\ln{e)\^ +CT , 

where c depends on k and a, but not e. Furthermore, the estimate \s £ — s§ | < Ck£ holds if Q 
is C 1 ' 1 (or convex in the case m = 1). 
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Remark 7.5. The operator S e introduced above is in fact the inverse of the Dirichlet-to- 
Neumann map associated with C £ . Note that by Corollaries 15.31 and 16. 2\ 



Sp — St 



Ce 



if n is C 



i,i 



o\\L 2 (an)^L 2 (dn) 



< 



C n e | ln(e) 1 2 +CT if fl is Lipschitz, 



for any a > 0. 



(7.2) 



8 Weighted potential estimates 



Let H £ {x) = (H](x), H™{x)) be defined by 

/- 

Jn dy k 



mix) 



{Tf(x,y)}h^y)dy, 



•1) 



where h = (h 1 , . . . , h m ) G L 2 (fi). It follows from [3] that || j| Z/2(Md) < C\\h\\ L 2 (n) . 
Proposition 8.1. The estimate 

1/2 



holds for any a > 1 . 



# e |U2(an) < C a < / 5(x)0 Q (5(a;)) dx 



Proof. Recall that 5(x) = dist(x, d£l) and 4> a (t) = {ln(| + e a )} a . Let g = (g 1 ,. 
L 2 (dQ) and u e = (u], . . . , u™), where 



4{y) 



d 



{Tf(x,y)}g a (x) da{x). 



It follows from [13, Theorem 4.3] that ||(«e)*IU 2 (dfi) — C|l IU 2 (dn)- Observe that 



H*(x)g a {x) da{x) 



no 



h\y)u^y)dy 



n 



< { ^ |/z( 2/ )| 2 < 5( 2 /)0 fi (5( 2 /))rf 2/ | 1/2 |^ I^Cg/)!^^^)^^^))}- 1 dy^ /2 



and that if a > 1, 



\u E (y)\ 2 {diy^MvW 1 dy < C / |(« £ )*| 2 da(y) < C\\g\\^ m) 

Jan 



Estimate (O follows from ([Q]) - (l8\Ij) by duality. 
Proposition 8.2. The estimate 



\H £ \\ L 2 {n) <C a \ ^ \h(x)\ 2 5(x)MS(x)) dx 



1/2 



12) 



•3) 



U) 
□ 

5.5) 



/io/ds /or any a > 1. 
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Proof. Let K C K\ be two compact subsets of Q such that dist(K, \ K\) > Co > 0. Since 
\W y T £ (x, y)\ < C\x — y\ 1 ~ d 1 we have 

\H e {x)\<C [ / |%)|dy foranya;Gir. 

J K\ \ x ~ y\ Jn\K! 

This implies that ||^ e ||L 2 (Ai is bounded by the right hand side of (j8.5p if a > 1. 

To estimate ||-ff e ||L 2 (Q\_ft:), it suffices to show that ||-ff e ||,L 2 (<9Q t ) is bounded uniformly in t 
by the right hand side of ( 18. 5 p for — c < t < 0, where Qt is defined in Section 2. This may 
be done by a duality argument, as in the proof of Proposition 18.11 Indeed, the argument 
reduces the problem to the following estimate 



\u £ , t (y)\ 2 {S(v)MS(y))} 'dy^C \g\ 2 da, (8.6) 

Jan t 



where 

<t(v)= f ^-{Tf(x,y)}g a (x)da(x) 



'an* fyjk 

Finally, the estimate (18. 6 p follows from the observation that ||w e ,t||L 2 (i<:) < C^r||fl'IU 2 (an t ) for 
compact K C Vt t , and that ||u e ,t||L3 (a n,) < C\\{u e j)*\\ L 2 {dQt) < C\\g\\ L 2 {dnt) for -c < t, s < 0. 
This completes the proof. □ 

Proposition 8.3. The estimate 

jjWH s {x)\ 2 5{x)<p a {5{x))dx\ <C a |y^|/i(x)| 2 5(x)0 a+2 (5(x))rfxj (8.7) 
holds for any a > 0. 

Proof. Using |V a; Vj / r £ (x, y)\ < C\x — y\~ d and a partition of unity, we may reduce the esti- 
mate (18. 7p to the case where Q = {(x', Xd) '■ x' G and Xd > ip{x')} is the region above a 
Lipschitz graph and S(x) is replaced by S(x) = \xd—ip(x')\. Since T e (x, y) = e 2 ~ d Ti(x/e, y/e), 
by a rescaling argument, we may further reduce the problem to the following weighted L 2 
inequality for a singular integral operator, 

f \T{f)\ 2 u x dx<C / \f\ 2 u 2 dx, (8.8) 
Jm d JR d 

where Ux(x) = 5(x)<p a (eS(x)), u^x) = 5(x)0 a+2 (e5(x)) and 

T(f)(x)= [ V x V y T 1 (x,y)f(y)dy. (8.9) 

JR d 

We point out that the constant C in (18. 8p should only depend on d, m, /i, A, r, a and 

HWHoo. 

To establish (18.81) . we first use the asymptotic estimates on V x V y ri(x, y) for \x — y\ < 1 
and | a; — y\ > 1 in [3] to obtain 

\T(f)(x)\ < C{\T*( 9l )(x)\ + \T*(g 2 )(x)\ + M(f)(x)}, (8.10) 
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where T*, T 2 * are L 2 bounded maximal singular integral operators with standard Calderon- 
Zygmund kernels, M(f) is the Hardy-Littlewood maximal function of / in M d , and \g\\, |^| 
are bounded pointwise by C\f\. Next we observe that u± is an weight in M. d . This allows 
us to use a classical result of R. Coifman and C. Fefferman [5] and (I8.10p to deduce that 



|T(/)|Wz<C / \M{f)\ 2 u l( ix. (8.1i; 
As a result, it remains only to show that 

|M(/)|V^ < C I \f\' 2 u 2 dx. (8.12) 



This is a two-weight norm inequality for the Hardy-Littlewood maximal operator, which has 
been studied extensively. In particular, E. Sawyer [21] was able to characterize all pairs of 
(wi,w 3 ) for which fl8A2|) holds. 

Finally, to prove (I8.12p . by a bi-Lipschitz transformation, we may assume that ip = 0. 
Consequently, it suffices to consider the case d — 1. This is because M(f) < Mi o M 2 o 
• ■ ■ o Mrf(/), where Mj denotes the Hardy-Littlewood maximal function in the Xi variable. 
Furthermore, by rescaling, we may assume e — 1. With ui(x) = \x\4> a {\x\) and uz(x) = 
|x|02+a(|^|) in h is not very hard to verify that (wi, U2) satisfies the necessary and sufficient 
condition in [21] for any a > 0. We omit the details. □ 

Proposition 8.4. The estimate 

\\M(HMLHan ) <C a yjh(x)\ 2 5(x)US(x))dx' ] j ' (8.13) 

holds for any a > 3. 

Proof. By the fundamental theorem of calculus and definition of the radial maximal operator, 
it is easy to see that for any Q G dfl, 



M{u){Q) <C y° {|Vu(A t (Q))| + \u(A t (Q)\}dt 

<C a y\\Vu(A t (Q))\ 2 +\u(A t (Q)\ 2 }\t\M\t\)dt^ ' 



(8.14) 



for any a > 1. This yields that 

I \M{u)\ 2 da < C a [ {\Vu(x)\ 2 + \u{x)\ 2 }5{x)(i) a {5{x))dx. (8.15) 
Jan Jn 

Letting u(x) = H £ (x) in (I8.15p . we obtain estimate (18.131) by Propositions 18. 2118731 □ 
Proposition 8.5. Let f = (f 1 , . . . , f m ) e L 2 (dVL) and u £ = (u l £) . . . , u™) be given by 

<(x)= [ 1 ^{Tf(x,y)}f(y)da(y). 

Then 

1/2 

|2 



Vu e {x)\ 2 8{x)dx\ <C\\f\\ L , m . (8.16) 
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Proof. By a partition of unity we may assume that Q = {(x', x d ) '■ x' G M. d 1 and Xd > ip{x')} 
is the region above a Lipschitz graph. By a rescaling argument we may further assume that 
e = 1. 

We first estimate the integral of |Viti(x)| 2 <5(x) on 

D = O + (0, . . . , 1) = {{x', x d ) : x d > V(a0 + !}■ 

By the asymptotic estimates of V x V ?; r 1 (x, y) for |x — y| > 1 in [3J, p. 906], we may deduce 
that if x G -D, 

|V, lW - W {x )\ < C f JM^M (8.17) 

Jan 1^ — 2/1 ' 

for some rj > 0, where is a finite sum of functions of form 

f 2 

e iA x ) / Q^-{ r f(^y)}^(y)da(y), 

with |ey(x)| < C and \g"\ < C\f\. Recall that Tq(x, y) is the matrix of fundamental solutions 
for the operator £ (with constant coefficients), for which the estimate (18. 16ft is well known 
[7j. It follows that 

|2 j„ ^ r> I I -f 1 2 



|W(a;)r<&(a:)dx < C / |/| 2 Ax. (8.18) 
u Jan 

Let /(x) denote the integral in the right hand side of (18 .171) . By the Cauchy inequality, 

i-v f \f(y)\ 2 da(y) 



|/(x)| 2 <C{8(x)}' 



an \x-y\ d+r > 



This gives J D \I(x)\ 2 5(x) dx < C\\f\\ 2 L2 ^ dn ^ and thus J D \V Ul (x)\ 2 5(x) dx < C\\f\\ 2 L2(m) . 
To handle V«i in \ D, we let 

A(r) = ((^'i il>(x')) : Ix'l < r), 

(8.19) 

T(r) = {(x,Xd) : \x'\ < r and ip(x) < Xd < ip(x) +Corj. 
We will show that if C\{u) = in the Lipschitz domain T{2) and (u)* G L 2 , then 

/ |Vu(x)| 2 |x d - tp{x')\dx < C I \u\ 2 da + C I \u\ 2 dx, (8.20) 

JT(1) JA{2) JT{2) 

which is bounded by C f A , 2 ) l( M )*| 2 ^ cr - By a simple covering argument one may deduce from 
(I8T20|) that 

/ \V Ul \ 2 5{x)dx < C [ |K)*| 2 Aj < C [ \f\ 2 d(j, (8.21) 
Jq\d Jan Jan 

where the last inequality was proved in [13]. 

Finally, to see (I8.20p . we use the square function estimate for L\ on T(r) for 3/2 < r < 2, 

/ |Vw(x)| 2 dist(x, <9T(r)) dx < C / \u\ 2 da, (8.22) 

JT(r) JdT(r) 
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to obtain 



/ \Vu{x)\ 2 \x d -^{x')\dx <C I \u\ 2 da + C [ \u\ 2 da. (8.23) 

JT(1) JA{2) JdT{r)\A(2) 

Estimate (18.201) follows by integrating both sides of (I8.23P in r G (3/2,2). We remark 
that under the condition A e A(/i, A,r), the square function estimate (I8.22p follows from 
the double layer potential representation obtained in [13] for solutions of the L 2 Dirichlet 
problem by a T(6)-theorem argument (see e.g. [HI pp. 9-11]). This completes the proof. □ 
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